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We investigate a one-dimensional S = 1/2 antiferromagnetic Heisenberg model coupled to 
quantum lattice vibration using a quantum Monte Carlo method. We study the ground-state 
lattice fluctuation where the system shows a characteristic structure factor. We also study the 
mass dependence of magnetic properties such as the magnetic susceptibility and the magnetic 
excitation spectrum. For heavy mass, the system shows the same behavior as the case of 
classical lattice vibration. On the other hand, for light mass, magnetic properties coincide with 
those of the static uniform chain. We investigate the physical mechanism of this behavior and 
propose the picture of quantum narrowing. 
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1. Introduction 

The spin-Peierls transition is one of the most fasci- 
nating features derived from cooperative phenomena of 
spin and lattice degrees of freedom due to the spin- 
phonon coupling in the one-dimensional antiferromag- 
netic Heisenberg model. The lattice shows a sponta- 
neous dimerization which accompanies the bond alterna- 
tion at low temperature, because a pair of spins tends to 
form a local singlet state on each strong bond to lower the 
magnetic energy which competes with the elastic energy 
of the lattice. An energy gap opens in the spin excita- 
tion spectrum. The magnetic susceptibility in all direc- 
tions exponentially drops to zero due to the spin gap as 
the temperature decreases. These phenomena have been 
observed in materials such as TTFCuS4C4(CF3)4^) and 
CuGeOg.^' 

Theoretical investigations have revealed ground state 
properties with an adiabatic approximation for the lat- 
tice. Usually, the exchange coupling is assumed to 
change in proportion to the lattice distortion amplitude 
e. Namely, the bond alternation is given in the form 
Ji — J[l + (— and 6 = ae. It has been shown 
by a bosonization method^' that the bond alternation 
brings about the magnetic energy gain proportional to 
j4/3 g _ ^12, which overcompensates the elastic en- 
ergy loss proportional to 5'^ . Thus the system dimerizes 
in the ground state for an arbitrarily small spin-phonon 
coupling. The adiabatic treatment for the lattice is valid 
when the mass of the magnetic ion is so heavy that the 
kinetic energy of the lattice vibration can be neglected. 
Put another way, the phonon energy ld is much small 
compared to the spin-Peierls gap and the exchange cou- 
pling. Recently we have studied finite temperature prop- 
erties of the adiabatic lattice case neglecting the kinetic 
energy term of the lattice by a quantum Monte Carlo 
method.^' We have observed how the bond alternation 



develops as the temperature decreases. 

On the other hand, quantum lattice fluctuation be- 
comes remarkable for the light mass case. It has been 
revealed that quantum lattice fluctuation disarranges the 
lattice dimerization in the ground state below a critical 
spin-phonon coupling. ^"^^ The quantum phase transi- 
tion between the dimerized phase and the uniform phase 
has been investigated with various methods such as an 
exact diagonalization method,^"' a density matrix renor- 
malization group method,^' ^^-^ a quantum Monte Carlo 
method^^"^'*-' and an analytical method.*' ^^-'^^'^^^ In ana- 
lytical investigations, the lattice degree of freedom is in- 
tegrated out for small w and the effect of the spin-phonon 
coupling is transposed to the dimerization and the geo- 
metrical frustration of the spin interaction in an effective 
spin Hamiltonian.^' Thermodynamic properties have 
been investigated by a quantum Monte Carlo method 
taking account of thermal fluctuation of the quantum 
phonon. ^"^^ For small lo, the lattice fluctuation can be 
regarded as the adiabatic motion and the lattice dimer- 
ization occurs at low temperature. There the magnetic 
susceptibility decays exponentially. In addition, it has 
been pointed out that thermal fluctuation of the lat- 
tice causes deviation of the magnetic susceptibility from 
that of the static uniform chain even at high temperature 
where the spin-Peierls correlation is not relevant. As w 
becomes larger, the magnetic susceptibility approaches 
that of the original spin chain of the static lattice, which 
suggests that spin and lattice degrees of freedom decou- 
ple and lattice fluctuation does not affect the spin state. 

In this paper, we pay attention to the mass depen- 
dence of the effect of quantum lattice fluctuation. In §2 
we explain a model and a numerical method. In §3 we 
study the mass dependence of the lattice configuration. 
For heavy mass, the lattice is uniform on the thermal av- 
erage at high temperature and dimerizes at low temper- 
ature. For light mass, the spin-Peierls correlation does 
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not appear even at low temperature. We find distinctive 
dependence of the structure factor on the mass. In §4 
we study the effect of quantum lattice fluctuation from 
a viewpoint of the world- line configuration of the lattice. 
In §5 we study the mass dependence of magnetic proper- 
ties such as the magnetic susceptibility and the magnetic 
excitation spectrum. For light mass, we find no effect of 
lattice fluctuation on magnetic properties, even though 
the lattice strongly fluctuates. In §6 we study the reason 
why the largely fluctuating lattice gives the same mag- 
netic behavior as that of the static uniform chain. It 
can be interpreted as a narrowing effect due to quantum 
lattice fluctuation. In §7 we summarize and discuss our 
results. 

2. Model and Method 

We investigate a one-dimensional S — 1/2 antiferro- 
magnetic Heisenberg model coupling with quantum lat- 
tice vibration, in which the lattice displacement produces 
the distortion of the exchange coupling between spins. 
The Hamiltonian is described by 

N 

H ^'^J[l + a{ui - Ui+i)] Si ■ Si+i 



(a) 



(b) 



i=l 
N r 



2m 2 



(1) 



where the distortion of the exchange coupling is assumed 
to be proportional to the lattice distortion. Here a is 
the spin-phonon coupling constant, m is the mass of the 
magnetic ion and k is the elastic constant. The periodic 
boundary condition is adopted for both spin and lattice 
degrees of freedom. In the simulation, we fix the total 
length of the lattice. We set the uniform exchange cou- 
pling J = 1 and take it as the energy unit. We consider 
a parameter set of a = 1 and k = 1 and investigate the 
behavior for various values of m. 

We use a quantum Monte Carlo method in order to 
investigate thermodynamic properties. We adopt the 
recipe introduced by Hirsch to deal with the spin-phonon 
coupled system.^*'' The partition function is expressed 
in the path-integral formula including both spin and lat- 
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Fig. 1. The schematic representation of the (1 + l)-dimcnsional 
lattice of AT = 6 and M = 2. 
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Fig. 2. The schematic representation of the change of the lattice 
configuration, (a) the local flip and (b) the global flip in which 
the lattice points denoted by solid circles are flipped together. 
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where M is the Trotter number and At — P/M with the 
inverse temperature (3. In Fig. 1 we show an example of 
the transformed (1 -I- l)-dimensional lattice. There are 
three layers in a unit for the imaginary-time axis, namely, 
two checkerboard layers for the spin-spin part and one 
layer for the phonon part. The unit of layers is heaped 
in M layers. 

We perform Monte Carlo simulations in the (1 -I- 1)- 
dimensional classical system. The spin configuration is 
updated by the loop algorithm^^' with the fixed lattice 
configuration. The loop algorithm is powerful to calcu- 
late physical quantities at considerable low temperature. 
Therefore it is useful to observe the spin-Peierls correla- 
tion at low temperature. In the loop algorithm, a graph 
element is assigned to each plaquette with a probability 
which depends on the value of the exchange coupling on 
the plaquette. In the present case, the value of the ex- 
change coupling varies along the imaginary-time axis due 
to quantum lattice fluctuation. Hence we can not take 
the limit of continuous imaginary time."'^^^ In this study, 
the lattice conflguration is updated by the Metropolis al- 
gorithm with the fixed spin conflguration. We perform 
a local flip and a global flip for the update of the lattice 
as shown in Fig. 2. The local flip updates the lattice 
state at a point in the (1 -I- l)-dimensional space, which 
causes a curved world line of the lattice configuration. 
For convenience of the calculation, the change of the lat- 
tice displacement is discretized by a small value Wunit- 
First, whether the lattice point moves right or left is 
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chosen with a probabihty 1/2. Next, whether the move 
takes place or not is chosen according to the Boltzmann 
weight. The local flip is efficient for the case where quan- 
tum lattice fluctuation is large. When quantum lattice 
fluctuation is small and the lattice tends to move adia- 
batically, we adopt the global flip of the lattice state. The 
global flip updates the lattice states by a parallel shift of 
all the points in a world line. We perform the spin update 
and the lattice update alternately. Therefore we obtain 
the thermal distribution of both spin and lattice degrees 
of freedom. Typically, the initial 10^ Monte Carlo steps 
(MCS) are discarded for thermalization, and the follow- 
ing 10^ MCS arc used to calculate physical quantities. 
The sampled data is divided into 10 bins, and the er- 
rorbar is estimated from the standard deviation for the 
data set of these 10 bins. 

For the simplicity of the notation, we use the bond 
distortion, 

Aj = a{ut - Ui+i), (3) 

as a variable instead of the lattice displacement Ui it- 
self. We provide a cutoff for the bond distortion |Aj| < 
Acutoff = 0.6, i.e., 0.4 < J, < 1.6, to conflne the ex- 
change coupling to be antiferromagnetic. As we have 
pointed out,^-* the variation of the bond is very large at 
flnite temperature if we use the Hamiltonian of eq. (1). 
Such large deviation is not realistic and we need some 
nonlinear suppression of the bond deviation. The cutoff 
plays the role of this suppression. 

3. Lattice Configuration 

The bond configuration characterizes the spin-Peicrls 
state, because the bond takes either a uniform configu- 
ration in the uniform state or an alternating configura- 
tion in the spin-Peierls state. We investigate the effect 
of both thermal fluctuation and quantum fluctuation on 
the bond configuration. 

In Figs. 3 we show the bond correlation function from 
the left edge bond, 



Cbond(l, j) = (AiA, ), 



(4) 



where (• • • ) denotes the thermal average. For a heavy 
mass m = 10000, the bond correlation function depicted 
in Fig. 3(a) agrees with our previous results without 
quantum fiuctuation of the bond.^^ At a high temper- 
ature T = 0.5, each bond fluctuates with no correla- 
tion and the bond takes a uniform configuration on the 
thermal average due to thermal fluctuation of the bond. 
At a very low temperature T = 0.02, the bond alterna- 
tion occurs with strong correlation extending over the 
chain. When the mass decreases, the effect of quantum 
lattice fluctuation becomes remarkable. As we show in 
Figs. 3(b) and 3(c), the alternating structure at low tem- 
perature disappears. We find no alternating pattern for 
a light mass m — 1, where the bond takes a uniform con- 
figuration because the bond alternation is disarranged by 
quantum fiuctuation of the bond. 

Thus the uniform bond configuration is realized due 
to thermal fluctuation at high temperature and due to 
quantum fluctuation at low temperature. The shape 
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Fig. 3. The bond correlation function along the real-space axis for 
(a) m = 10000, (b) m = 100 and (c) m = 1. The left figures are 
the data at T = 0.5 for the system of iV = 64 and M = 12. The 
right figures axe the data at T = 0.02 for the system of AT = 64 
and M = 384. 



of the bond correlation function is simply flat for both 

cases. Therefore the bond correlation function is insuf- 
flcient to characterize the difference of the nature of the 
bond fluctuation in both cases. For the purpose of clari- 
fying the distinction between the classical bond fluctua- 
tion and the quantum one, we investigate the bond struc- 
ture factor. 



where 



.5bond(«) = (|A,|2), 



1 ^ 



(5) 



(6) 



In Fig. 4(a) we show the bond structure factor for m = 
10000. At T = 0.5, S'bond is flat as a result of the mix- 
ture of various modes of the motion, which means that 
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Fig. 4. The bond structure factor for (a) m = 10000, (b) m = 100 
and (c) m = 1. The left figures are the data at T = 0.5 for 
the system of AT = 64 and M = 12. The right figures are the 
data at T = 0.02 for the system of AT = 64 and M = 384. The 
open circle in (c) denotes the data at T = 2.0 for the system of 
AT = 64 and JM" = 4. 



each bond fluctuates independently. At T = 0.02, the 
lattice forms the bond-alternate structure and /Sbond has 
a Bragg peak at g = tt. Except for q = n, ^bond is very 
small and the fluctuation is independent of the position 
of the bond. These structures are consistent with the 
results of the classical lattice system. 

On the other hand, when the mass decreases, the bond 
structure factor shows characteristic dependence on the 
temperature. In Fig. 4(b) we show the bond structure 
factor for m = 100. 5bond remains flat at T = 0.5. At 
T = 0.02, comparing with the data for m = 10000, the 
amplitude of the Bragg peak becomes small and the fluc- 
tuation of various modes grows in addition to g = tt. This 
structure factor with the mixture of various modes cor- 
responds to the correlation function shown in the right 
of Fig. 3(b), where the alternate structure is hardly de- 



tected. 

When the mass decreases furthermore, m = 1, S'bond 
has a sinusoidal shape at both T = 0.5 and T = 0.02 
as shown in Fig. 4(c). The sinusoidal structure factor is 
interpreted as follows. For light mass, the second term 
of eq. (1) is much larger than the first one. Namely, the 
lattice moves almost independently of the spin configu- 
ration. Therefore, in the following, we investigate the 
lattice part of eq. (1), 



N 

.^phonon ~ ^ ^ 



(7) 



The operators of the lattice displacement Uj and the mo- 
mentum pj are expressed by the creation operator aj and 
the annihilation operator Uq of the phonon. 



Pj = 



where 



2 A 2 ■ 2 9 
1 9 



Then eq. (7) is transformed into 

■ffphonon = '^^q + 2^ " 



(8) 
(9) 

(10) 

(11) 



The Fourier component of the bond distortion is ex- 
pressed by 



Aq = 2iai 



1 



2Nmu>„ 



(aj -I- a-q)e 



-'«/2sin|. 



(12) 



At r = 0, all the nomal modes are in the ground state, 
i.e., aq\0) = 0, and ^bond is given by 



S'bond(9) 



(|A,|^) 

2a2 



a q 
sm - . 

Nmuj 2 



(13) 



Namely, the zero-point motion gives the sinusoidal struc- 
ture factor. This result agrees with the sinusoidal shape 
of our Monte Carlo results. But we find that the am- 
plitude in Fig. 4(c) (= 0.0021) is smaller than that of 
eq. (13) (= 0.0156). This discrepancy is caused by the 
cutoff for the bond distortion. The bond distortion is 
suppressed due to the cutoff and the structure factor be- 
comes small. In Monte Carlo simulations of the lattice 
part without the spin part, it is confirmed that the am- 
plitude agrees with that of eq. (13) in the case where the 
cutoff is large enough, and the amplitude agrees with 
that of Fig. 4(c) in the case where the cutoff is the same 
as that we provide in the spin-Peierls model. Thus we 
conclude that the lattice fluctuation for light mass is 
characterized by the zero-point quantum fluctuation of 
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the lattice. It should be mentioned that when the tem- 
perature is larger than w = ^Jk/m, ^bond is influenced 
by the excited state of the phonon. In the present case, 
we actually find that S'bond is influenced by excited states 
atT = 2.0 (>a; = 1.0) and changes from the sinusoidal 
shape to a flat shape as shown by the open circle in the 
left of Fig. 4(c). 

4. Quantum Lattice Fluctuation 

In the thermal average, the lattice takes either the uni- 
form configuration for light mass or the dimerized config- 
uration for heavy mass due to the spin-phonon coupling. 
On the other hand, in a snapshot in the Monte Carlo 
simulation, the world line of the lattice shows various 
configurations due to quantum lattice fluctuation. For 
heavy mass, the motion of the lattice is regarded to be 
adiabatic and the world-line conflguration of the lattice 
is straight along the imaginary- time axis. The straight 
world line parallels with each other and forms the dimer- 
ized configuration. As the mass decreases, quantum lat- 
tice fiuctuation becomes relevant and the world line of 
the lattice begins to fluctuate along the imaginary-time 
axis. The deviation of the world line may suppress the 
bond- alternating order. Here we introduce the order pa- 
rameter for the bond-alternating order. 



(a) 0.1 



A2 

sg 



(14) 



In addition, we explore how the bond alternation is dis- 
arranged due to quantum lattice fluctuation by investi- 
gating the following quantity. 




N 



^(-l)^e^-A, 



-Ht 



where e^'^Aj e~^'' 
imaginary time r. 

dencc of AL and 



Agg agrees 



of A^g 

with 



A2 , 
AL 



(15) 

represents the value of Aj at the 
In Figs. 5 we show the size depen- 
at T = 0.02. For to = 10000, 
which indicates that the bond- 
alternating order develops in both the real-space axis and 
the imaginary-time axis. In fact, there the world line is 
straight along the imaginary-time axis. When the mass 
decreases, the absolute values of the bond-alternating 
correlations decrease as shown in Figs. 5(b) and 5(c). In 
particular, A^g^. is much reduced comparing with A^g for 
flnite size systems. This is due to the deviation of the 
world-line configuration along the imaginary-time axis. 
For TO = 100, however, we find that A^g and A^g^ con- 
verge to almost the same finite value in the thermody- 
namic limit. The amount of the bond alternation for 
TO = 100 is smaller than that for to = 10000. For to = 1, 
the bond configuration is uniform on the thermal aver- 
age in the thermodynamic limit. There A^g shows 1/N 
dependence which indicates a short range order of the 
bond alternation, while A^g^ is almost zero because the 
correlation length along the imaginary-time axis is much 
shorter than (3 as we will see in §6. 
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Fig. 5. The size dependence of the order parameters for the bond 
alternation A^g (solid circle) and A^gj (open circle) for (a) m = 
10000, (b) m = 100 and (c) m = 1 at T = 0.02. M = 384. 



5. Mass Dependence of Magnetic Properties 

The spin-plionon coupling causes deviation of mag- 
netic properties from those of the uniform lattice sys- 
tem. First we consider the effect of thermal fiuctuation 
of the lattice at high temperature where the spin-Peierls 
correlation is not relevant. In Fig. 6 we show the tem- 
perature dependence of the magnetic susceptibility for 
a high temperature region. We have calculated for the 
adiabatic lattice system. For comparison, the data for 
the adiabatic lattice system (dotted line) and the uni- 
form lattice system (solid line) are shown in the figure. 
For to = 10000, the magnetic susceptibility agrees with 
that of the adiabatic lattice system. As has been pointed 
out in numerical works, ^' ^^-^ the magnetic suscepti- 
bility deviates from that of the uniform lattice system 
even at high temperatures due to thermal fluctuation of 
the lattice. As the mass decreases, the magnetic sus- 
ceptibility approaches that of the uniform lattice system 
monotonously. For to = 1, the magnetic susceptibility 
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Fig. 6. The temperature dependence of the magnetic susceptibil- 
ity. The data are calculated for the system of iV = 64 and 
extrapolated with respect to the Trotter number. The dotted 

line denotes the data for the adiabatic lattice system of A'' = 64 
and the solid line denotes the data for the uniform lattice system 
of AT = 64. 
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Fig. 7. The magnetic elementary excitation spectrum in the 
ground state. The data are calculated at T = 0.02 for the sys- 
tem of iV = 64 and M = 384. The sohd hne denotes the data 
for the uniform lattice system oi N = 64. 
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agrees with that of the uniform lattice system weh, even 
though the lattice strongly fluctuates. This fact has been 
seen in previous works for models with the second quan- 
tization of the lattice. 

Now we consider the effect of quantum fluctuation of 
the lattice at low temperature. In order to grasp the 
magnetic behavior from a microsc;opic; point of view, we 
investigate the magnetic excitation spectrum. It can be 
extracted from the imaginary-time correlation function 
of the spin,^^' 

SsMi, r) = {e"^SIe-"^Sl^), (16) 

where 

1 ^ 

The dispersion relation E{q) of the elementary excitation 
is obtained by 

E{q) = -^\n[S,Mq,T)] (18) 
for large r. In Fig. 7 we show the magnetic elementary 



Fig. 8. The imaginary-time correlation functions of the spin and 
the bond for <? = tt at T = 0.02. (a) m = 10000, (b) m = 100 
and (c) m = 1. The solid circle and the open circle denote the 
data in terms of the spin and the bond, respectively. The data 
axe calculated for the system of = 64 and M = 384. 



excitation spectrum at T = 0.02 for several values of the 
mass. The solid line denotes the data for the uniform lat- 
tice system of = 64, where we find a small energy gap 
due to the finite size effect. For m = 10000, the lattice 
motion is adiabatic and the bond alternation occurs at 
low temperature. Accordingly, we find an energy gap at 
q = TT clearly. On the other hand, we find that the energy 
gap becomes small drastically as the mass decreases. For 
m = 1, E{q) coincides with that of the uniform lattice 
system. Thus we find that not only the magnetic suscep- 
tibility but also the microscopic magnetic characteristic 
agrees with those of the uniform lattice system. There- 
fore we conclude that this agreement is intrinsic. In the 
next section, we study the physical background of this 
agreement. 
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6. Quantum Narrowing Effect 

In this section we consider the reason why the expecta- 
tion vahic of magnetic properties of the light mass system 
coincides with that of the uniform lattice system. As we 
mentioned above, the motion of the lattice is almost in- 
dependent of the spin configuration because the mass is 
light. On the other hand, the motion of the spin may 
be affected by the bond fluctuation. Generally speak- 
ing, however, the effect of fast fluctuation is to be ne- 
glected, which is known as the motional narrowing effect. 
Namely, when the fluctuation is much fast compared to 
the time scale which we would like to observe, the fluc- 
tuation is leveled off and does not affect other degrees of 
freedom. In the present case, what comes into question 
is the degree of quantum fluctuations of the spin and the 
lattice. The correlation along the imaginary-time axis 
gives information on the time scale of quantum fluctua- 
tion. We investigate Sspin{q,T) defined in eq. (16) and 
the imaginary-time correlation function of the bond, 



5bo„d(a,T) = (e«"Age-^^A_,). 



(19) 



In Figs. 8 we compare both imaginary- time correlation 
functions for g = tt at T = 0.02. For m = 10000, 
Sspin{q,T) decays along the imaginary-time axis quickly 
and <S'bond(5, t) keeps a constant value. The spin fluctu- 
ation is much faster than the bond fluctuation. Namely, 
the time scale of the spin fluctuation Tgpin is much smaller 
than that of the bond fluctuation Tbond ("^spin -C Tbond)- 
In this case, the lattice vibration can be regarded to be 
adiabatic. On the other hand, as the mass decreases, 
the inequality relation between the two time scales is re- 
versed. For m = 1, S\)ond{<l,T) decays quite rapidly. This 
value of Tbond is much shorter than /3 = (0.02)~^ = 50 
as we expected by the observation of A^g^ in §4. Here 
the relaxation of Shond{q,T) is much faster than that of 
5'spin(g,r), i.e., Tspm > TTjond- This means that the lat- 
tice fluctuates much faster than the time scale of the 
motion of the spin. Thus we conclude that a sort of nar- 
rowing effect happens and the lattice is effectively fixed 
to be uniform. We call this phenomenon "quantum nar- 
rowing effect" . 

7. Summary and Discussion 

In this paper we investigated a one-dimensional S = 
1/2 antiferromagnetic Heisenberg model coupled to 
quantum lattice vibration using a quantum Monte Carlo 
method. We studied the mass dependence of various 
physical quantities. For heavy mass, quantum lattice 
vibration is regarded as the adiabatic motion and the 
system dimerizes at low temperature. When the mass 
decreases, the bond-alternate correlation is reduced due 
to the deviation of the world-line configuration of the lat- 
tice. For light mass, the spin-Peierls correlation does not 
grow even at low temperature, and the ground-state lat- 
tice fluctuation is equivalent to the zero-point vibration 
of the phonon. The lattice fluctuation does not affect 
the magnetic behavior. We found the separation of time 



scales of quantum fluctuations of spin (rspin) and lattice 
(Tbond) degrees of freedom. The time scales change with 

the mass. The motion of the lattice can be regarded 

to be adiabatic for heavy mass, where Tspin ^ Tbond- 
On the other hand, the motion of the lattice is smeared 
out due to the narrowing effect for light mass, where 

Tspin Tbond- 

For intermediate mass, ni = 100 [Fig. 8(b)], 
imaginary-time correlation functions of the spin and the 
bond decay with almost the same decay time. In such 
cases we expect that the spin fluctuation and the bond 
fluctuation resonate with each other. The resonance of 
quantum fluctuations would cause a multiple quantum 
state of spin and lattice degrees of freedom, which is an 
interesting problem in the future. 
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